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Abstract 

We investigate the influence of an external magnetic field on spin phase re- 
laxation of single electrons in semiconductor quantum dots induced by the 
hyperfine interaction. The basic decay mechanism is attributed to the dis- 
persion of local effective nuclear fields over the ensemble of quantum dots. 
The characteristics of electron spin dephasing is analyzed by taking an aver- 
age over the nuclear spin distribution. We find that the dephasing rate can 
be estimated as a spin precession frequency caused primarily by the mean 
value of the local nuclear magnetic field. Furthermore, it is shown that the 
hyperfine interaction does not fully depolarize electron spin. The loss of ini- 
tial spin polarization during the dephasing process depends strongly on the 
external magnetic field, leading to the possibility of effective suppression of 
this mechanism. 
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The spin state of an electron confined in a semiconductor quantum dot (QD) is considered 
one of the most promising candidates for reahzing the basic building block (i.e., qubit) of 
a quantum information system. Since the fundamental concept of this new paradigm relies 
on quantum mechanical entanglement of qubits, it is quite crucial to control spin relaxation 
processes that destroy the coherence of spin quantum state. So far, most of the attention 
has been devoted to the relaxation processes that result in irreversible loss of wave function 
phase due to spin-phonon interaction caused by spin-orbital coupling in solids or hyperfine 
interaction (HFI) in crystals with non-zero nuclei spin moments (see Refs. |l] and|^ as well as 
the references therein). The common feature of these spin-lattice mechanisms is that their 
relaxation rates are very small in QDs at low temperatures. 

On the other hand, the HFI can be considered as a source of a local magnetic field H up 
acting on the electron spin that does not disappear at low (or even zero) temperature.! 
This particularity makes the HFI a potentially dominant mehchanism at sufficiently low 
temperatures. In typical QDs, a sum of the contributions from a great number of nuclei 
spins forms this field. Thus, the strength and the direction of the H hf are the random 
variables, which vary from one QD to the next. Obviously, this dispersion can be damaging 
to quantum computation since electron spin precession occurs with a random phase and 
frequency. Nevertheless, it appears from a qualitative speculation that the role of H hf 
dispersion diminishes progressively with an increasing strength of the homogeneous external 
magnetic field B applied to the array of QDs. 

In this paper, we provide a quantitative analysis of electron spin evolution under the 
presence of an external magnetic field B as well as the local hyperfine field H hf- Note, 
that a theory of electron spin relaxation caused by HFI in a QD was recently presented in Ref. 
^. The main interest of Khaetskii et al% lies on the electron spin decoherence process inside 
a single QD when the external magnetic field is zero. It also contains a brief discussion of 
spin dephasing time. On the other hand, the present study concentrates on the important 
process of electron spin dephasing induced by HFI in an ensemble of QDs and explicitly 
considers the effect of external magnetic fields. 
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The Hamiltonian of electron spin S* in a QD containing nuclear spins Ij {j = 1, A^) 
takes the form {h = 1) 

N N ^ 

H = UeS^+ Un^Ijz + ^Aj S ■ I j, (1) 

i=i i=i 

where uJe and ujn are electron and nuclear spin splitting in a magnetic field directed along 
the Z axis, Aj is a constant of HFI with the j-th nuclear spin. The Hamiltonian of Eq. 
(|^) is isomorphic to that, which was introduced to describe the bound magnetic polaron in 
diluted magnetic semiconductors. Most particularities of optical spectroscopy as well as the 
thermodynamics of bound magnetic polaron were successfully described by a model where 
differences in constants of spin-spin interaction were ignored.lii By analogy with this model, 
we describe the dynamics of electron spin by Hamiltonian with a single effective constant 
A of HFI. This allows one to express the Hamiltonian in terms of the total nuclear spin 
moment F= H^=i I- 

n = UeSz + OJnF, + A 5^ ■ F . (2) 

Even with such a reduced Hamiltonian, the corresponding equations for the motions of 
electron and nuclear spins are quite complex, 
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^^[S)= B,x{S) + A(^F xSy, (3) 

j^{F)^B.x{F)^A{sxF). (4) 

where B e = {0, 0, oOe} and i? „ = {0, 0, oOn} are in units of energy. Assuming that the fiip-fiop 
processes are unimportant for the problem we consider, the ^'-function of the Hamiltonian 
given in Eq. can be factorized with respect to electron and nuclear spins. This means 



that ( F X S 



F)x(S 



and Eqs. and (||) present the closed system of equations 
with respect to {So) and (Fa) (a = x, y, z). 

Despite the seeming similarity of Eqs. (|^) and (^, there is a significant quantitative 
difference in the effective fields of the nuclei acting on the electron spin AF and of the 



electron ^Aj {S = 1/2) acting on the nuclei since a large number >> 1 of nuclear spins 
are involved in the QD. The latter property with regard to the inequality Ue » uJn means 
a large difference in the typical precession periods for electron and nuclei spins. Thus, for 
a moment, we can consider the time-dependence of ^ F ^ that is given in a parametric 
representation. The corresponding solution of Eq. with respect to X(t) = {Sx), y(t) = 
{Sy), Z{t) = (Sz) and initial conditions X(0) = F(0) = 0, Z{0) = 1/2 reads 

AF AF 
X{t) = -^{AF, +uj,){l - cosQt) + ^ sin fit; (5) 

AF AF 
Y{t) = i^{AF, + ^e)(l - cos fit) - ^ sin fit; (6) 

1 A^iF'^-F'^) 
m = - - - cosfit). (7) 

The frequency of electron spin precession 



fi = ^A^F^ + 2AF,uj^ + uj^ (8) 

has the simple physical meaning of the Zeeman frequency of electron in the total field 
H e = B e + AF composed of external and internal nuclear fields. 

Actually, the components of the total nuclear spin moment Fa(t), (a = x, y, z) is affected 
by the external magnetic field and the effective field of an electron aI s\ which is 



oscillating with a high frequency fi near some mean value Ai^S j. This mean value also 
changes slowly with a typical nuclear frequency fi„^ B n + A/2. By taking a small interval 
in time At (<< fi^^) to be longer than the period of electron beats (i.e.. At >> fi^^), one 
can average over the rapidly oscillating electron contribution in Eq. and arrive at a 
system of equations with respect to F^ = Fa{t): 
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-F, = -uJnFy + AY{t)F, - AZ{t)Fy- (9) 



-Fy = oOnF, + AZ{t)F, - AX{t)F,; (10) 



^F, = AX{t)Fy - AY{t)F,. (11) 

Here the bars mean the procedure of setting sinf2t and cos fit to zero when Eqs. (|])-(^ are 
substituted into Eqs. (|^)-([TT|). 

Clearly, Eqs. (p|)-([TT|) are non-linear. Fortunately, we obtain identically zero for the time- 
derivative Fz= that provides conservation of the nuclei spin projection on the external 
magnetic field during the dephasing process described by Eq. (|^). Thus, the dispersion 
of F = F and Fz, caused by the thermal fluctuation in an ensemble of nuclear spins 
associated with QDs, results in the dephasing of electron spin. 

Mathematically, evolution of the spin polarization in an ensemble of QDs can be reduced 
to the problem of averaging Eq. (|^) over the distribution functions for total nuclear spin 
[P{F)] and for its ^-projection = F^ [P^ (//)]: 

{Z{t))^ = J J Z{t)Pzifi)P{F)d^idF (12) 

The corresponding distribution functions were found in Ref. |^ (see also Ref. 0). In the 
case of non-saturated nuclear spin polarization, Eq. (0) can be well approximated by the 
expression 

{Z{t))^ = C dF j_^d^^FeM-— " ^)^(^), (13) 

where a = |/(/ + 1)A^ and C is a constant. The averaging in Eq. (|T3|) can be performed 
numerically for different experimental situations. Some examples of these calculations are 
presented in Fig. 1 for the case of relatively high temperature T ^ AauJn/we. 

Two important conclusions follow immediately from Eqs. (^ and (^) (see also Fig. 
1). First, electron spin dephasing does not fully relax its own initial spin polarization 
Zq = Z{0). The minimal rest Z^o = limj^oo {Z{t))rp of Zq corresponds to zero magnetic 
field: Zoo = |-^o- In a strong enough magnetic field {B oo), dephasing will be fully 
suppressed: Z^o Zq. Intermediate cases can be traced with high enough accuracy if we 
substitute high-temperature approximations for the average values of yU and yU^ in Eq. (|13[): 



~ I (F^) ^p- = ^ and (/U^)^ ~ ^ (F^) = |. A more detailed approximation is of no 
interest since it corresponds to saturation of nuclear spin polarization along electron spin 
polarization that results in insignificant dephasing. The final result reads 
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- a{\A- + Au^-^) + uf'- ^^^^ 

The second conclusion is that dispersion of Zeeman frequency Vt controls the electron 
dephasing rate. Hence, instead of calculating {Z{t))j, from Eq. ([13|) , we can find the dis- 
tribution function for Vt] then, the dephasing rate can be estimated as the width of this 
distribution. In other words, the dephasing rate can be obtained from the width of electron 
spin resonance lineshape 

g{uj) = j j 6{uj- n)P,{i2)P{F)dfidF. (15) 

Using the definitions in Eqs. ([T3|) and (|]) as well as the new integrand variables x = 
F/y/a, y = fJ^/y/o-, Eq. ( [T5|) can be rewritten as 

g{uj) = C'ooe^'^ "^o sinh {2LJeTo + \/o'uJn/T)ujTo . (16) 

where Tq = l/A^o". The dephasing rate t^^ is expected to be the half width of the 1/e 
decay in the maximal intensity of g{uj). Equation ([16|) predicts a weak dependence of t^^ 
on the parameter 2uJeTQ + ^/o^uJn/T, which is proportional to the external magnetic field: 

r^^ = (p{2LJeTo + A/au;„/T)ro'\ (17) 

where the function if{x) falls into the region between (0.69,1) (see Fig. 2). Thus, by the 



order of magnitude, r^^ c:^ = A^^I{I + 1)A'". 

At the first glance, it would seem that the efficiency of the phase relaxation increases 
with the QD volume Vq as y/N = VrijVo, where rij is the concentration of isotopes with the 
nuclear spin /. In actuality, however, we should take into account that the constant of the 
contact interaction A is proportional to the electron spin density at a nuclear site.i As a 
result, the rate t^^ must reveal an inverse dependence on ^/Vo. 



Let us estimate the HFI constant A in terms of ENDOR experimental data in Si:P.li 
According to the definition, 



A = -^geQn^^BfJ'n ^ (^n) 



(18) 



Here Qe and Qn are electron and nuclear g-factors, fiB and fin = {fT^e/^p) f^B the Bohr and 



nuclear magnetons, rrie and rrip the electron and proton masses, \E' ("^"j electron 
envelope function density at a nuclear lattice site T*„, and the parameter rj refiects the 
enhancement of electron density due to the Coulomb singularity at the nuclear core [see, for 
details, Ref. ^ and Eq. (2.78) therein]. In the case of Si^^, this parameter is r/ = 186. 

To find the parameter A for QD electrons in natural Si, we can use the Feher's exper- 
imental datai obtained for a shallow donor electron in the presence of Si^^ separated by 
[400]a from the donor centrum (4a = 5.43 x 10"^ cm): A^oof^T^^ = 7.72 MHz. At the same 
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time, the donor electron density ^ (^"^'400 j calculated at this location is 0.4 x 10^^ cm~^. 
Then, by utilizing the functional form of Eq. (|18[) that is proportional to the square 
of the electron density, the HFI constant A at a nuclear site can be expressed as 



r 400 ) • If a uniform envelope function is assumed inside the 
QD, \E' ("^n) ~ l/v^K)5 can find A/27Th = 1.93 kHz with the QD diameter and thick- 
ness of 500 A and 50 A, respectively. 

Then, in terms of the uniform \E'-function approximation, we can estimate the number of 
isotopes Si^^ in a QD: = xVq/VIo = 2.29 x 10"^, where x = 0.0467 is abundance of isotopes 
Si^^ and Qq = 20 the volume per one Si atom. Thus, one can find a = |/(/ + 1)A^ = 
1.14 X 10^ and ro=0.54 x 10"^ s. 

Let us also consider the final loss of electron polarization after dephasing reaches a 
saturation at t ^ To. According to Eq. (|14D, the relative loss of initial electron polarization 
Zq is 
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(19) 



Zo 3/2 + (l + £)(a;ero)2 

where e = (aA/T) (un/uje)- The value of aA which is independent of electron radius local- 
ization and abundance of isotopes Si^^ estimates to be approximately 0.025 K. Thus, with 



the exception of extremely low temperatures, we expect the parameter e to be small. 

We can see that Eq. (p!9| ) looks similar to that of Hanle effecti but with the opposite 
meaning: increasing magnetic field leads to enhancement (not suppression) of spin polar- 
ization. It reflects the competition of random field of nuclei, Hhf = ^/gef^BTo (that is on 
the order of 0.1-1 G in the case of Si QDs), and the external magnetic field B. Accordingly, 
dephasing is not important if i? ^ Hhf- It is interesting to note that Hhf can reach 
hundred G in QDs composed of III-V semiconductors.! 

In conclusion, we have considered the influence of an external magnetic field B on the 
electron phase spin relaxation in QDs. The dispersion of nuclear effective local fields H hf 
over the ensemble of the QDs results in electron spin dephasing. The dephasing rate can 
be estimated as a spin precession frequency caused by field H hf, so it depends weakly 
on an external magnetic field B provided nuclear spin polarization (or magnetization) is 
not saturated. On the other hand, reduction of electron spin polarization along B during 
the dephasing process does not lead to full depolarization. If the magnetic field is suitably 
strong, B ^ Hhf, one can expect an insignificant electron spin depolarization caused by 
the dephasing process considered. 

The result obtained in this study is directly relevant to quantum computing with a large 
number of qubits (i.e., a large number of QDs containing single electrons). Clearly, random 
variation in the spin precession frequency and phase across the qubits is deterimental; for 
example, it can break/negate the initial state preparation in a relatively short time frame 
of approximately tq. Special methods must be applied to compensate this distortion. Our 
work suggests an effective way to suppress the dispersion of random local fields, making QDs 
a more attractable candidate for quantum computing. 

The work performed at North Carolina State University was supported, in part, by the 
Office of Naval Research and the Defense Advanced Research Projects Agency. 
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FIGURES 

FIG. 1. Spin dephasing of QD electrons calculated for three values of magnetic fields {oOgTo =0, 
1, and 2). The temperature is assumed to be relatively high: i.e., T 3> AauJn/i^e (non-saturated 
nuclear spin polarization). 



FIG. 2. Calculated value of function 4'{x). 
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